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Abstract. The repulsive Coulomb force poses severe challenges when describing (d, p) 
reactions for highly charged nuclei as a three-body problem. Casting Faddeev-AGS equa¬ 
tions in a Coulomb basis avoids introducing screening of the Coulomb force. However, 
momentum space partial-wave f-matrix elements need to be evaluated in this basis. When 
those t-matrices are separable, the evaluation requires the folding of a form factor, de¬ 
pending on one momentum variable, with a momentum space partial-wave Coulomb 
function, which has a singular behavior at the external momentum q. We developed 
an improved regularization scheme to calculate Coulomb distorted form factors as the 
integral over the Coulomb function and complex nuclear form factors. 


1 Introduction 

Direct reactions involving rare isotopes offer an important tool for understanding the structure of such 
nuclei. Deuteron induced reactions are attractive from the experimental perspective, since deuterated 
targets are readily available. Theoretically they are attractive, since the scattering problem can be 
reduced to an effective three-body problem, which can be solved exactly using Faddeev techniques. 
The momentum space Faddeev equations (here in the Alt-Grassberger-Sandhas form) have success¬ 
fully been solved for {d, p) reactions involving light nuclei IH]. However, the screening technique em¬ 
ployed to handle the Coulomb force encounters technical difhculties, when applied to {d, p) reactions 
with heavier nuclei An alternative to the screening procedure is a solution of the Faddeev-AGS 
equations in the Coulomb basis. This was suggested and carried out in Ref. I^l using real two-body 
transition operators in separable form, neglecting spin degrees of freedom. 

Casting Faddeev-AGS equations in the Coulomb basis requires the evaluation of momentum-space 
partial-wave matrix elements of two-body transition operators. 



( 1 ) 


^e-mail: to@vsl.name 
^e-mail: elster@ohio.edu 



EPJ Web of Conferences 


If ?;(£’) is separable, one needs to consider 

t^{q, q\ (2) 

zy 


In this case only the folding of a form factor ui{p), depending on one momentum variable, with 
partial-wave Coulomb function for which we developed the numerical procedures in Ref. 0, 

is required: 

= J and uf/q'Y = J (3) 

Both integrals have an oscillatory singularity in the point p - q, where q is the external momentum. 
The partial-wave Coulomb functions are given by 




Six, n) 

(X + i0)‘+'’' 


!B(x,q,l,qy 
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X = p-q- (4) 


To evaluate u^J^q), one needs to regularize the integral in Eq. (I3]i. In Ref. fl this regularization is 
performed by using a Gel’fand-Shilov (see Ref. |01) technique for real form factors. When working 
with form factors describing complex nucleon-nucleus potentials, the regularization scheme must be 
constructed for complex form factors. 


2 Regularization Procedure 


Here we present a regularization scheme for calculating Coulomb distorted form factors as inte¬ 
grals over a Coulomb function and a complex form factor, which differs from the scheme presented 
in Ref. ||2l. The key point of the work of Ref. ||3l is that the integral is regularized only in a tiny 
neighborhood (A ~ 10“^ fm ') around q. That allows that the regularized singular part of the folding 
integral can be reduced to an analytic expression. However, the non-regularized parts of the inte¬ 
gral must be calculated quite close to the value of q, which requires a lot of integration mesh points 
independent of the numerical integration quadrature used. For example, we used thousands of Gauss- 
Legendre quadrature points to obtain results converged to 3 significant figures. 

To avoid using an excessive number of integration points, we developed here the new version 
of the regularization scheme, where we apply the regularization in a small, but finite neighborhood 
(A ~ 10“^ fm ') around q. Now the non-regularized integrals are converging much faster to an even 
better accuracy (see Fig. [T]i, while the regularized part is still mostly dominated by the analytical 
terms. Despite the fact that we need to calculate the regularized integral, we do not need to calculate it 
as accurately, as we had to calculate the non-regularized parts in Ref. ist]. Thus, the required amount 
of CPU time is greatly reduced. 

Similar to Ref. ISD, we start by splitting the folding integral into regularized and non-regularized parts. 


jf r /'A _ rt'-A , rt'+A , TOO 




We first focus on the regularized part which we rewrite by substituting Eq. (|4]i into Eq. (|5ll 
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Figure 1. The error plot of the non-regularized 
integrals /* with the same lower limit a = 0.4 fm“‘, 
but different upper limits: [solid 
black] b = 0.499 fm^' (i.e. A = 1 • 10^^ fm^*), 
[dashed red] b = 0.499999 fm-‘ (A = 1 • lO^*^ fm-‘). 
The parameters are: q = 0.5 fm“*, I = 0,q = 5.3. 
These values of q and q correspond to the p + ^°*Pb 
system with F^.m. = 6 MeV. See Sects, [^andj^for 
discussion and definition of the relative error cr*. 


where 


(f>+{x,q',i,T]) = ui,ij{q' +x)iq' +xf^(x^q'^lq)^ 

(p-(x,q',i,q) = ui,y(q' + x)(q' + xf^(x^q'■ 

To apply the regularization scheme, we split the integrals J± and take the limit +i0, 
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While taking the limit, we used the following property. 

By applying the Gel’fand-Shilov regularization technique from Ref. |01, we obtain 


(7) 
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where 
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( 12 ) 

(13) 


To calculate (p'^iO, q', l,n), we also need duiy{p)ldp. Taking the limit A —> +0, one arrives at the 
expressions from Ref. 

By using the expressions from Eq. (fTTT i and Eq. (|6]l, we can compute the regularized part for 
any given finite A. Then, using Eq. (|5]l we evaluate uf (q')^, while uf (q) is given by 

‘^y hZ 
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To calculate (p'JiQ, q', I, q), the expression for dSjdx is required. By using Eq. (B9) from Ref. ||3l, 
since p - q' + x, Eq. (15.5.1) from Ref. ||3l, and by utilizing the property of the hypergeometric 
function 2 Fi {...; 0) = 1, we obtain 


dSix, q, I, q) 


dx 




= e-“^T(l+/77)(-3 + 


(15) 


Here we substituted q' by q to simplify the expression. 


3 Results, Discussions and Outlook 

To prevent the unreasonably large computational times for calculating Coulomb distorted form factors, 
we improved the regularization scheme of Ref. Ist]. This new scheme converges much faster, as the 
number of integration points increases. To demonstrate it, we present the error plot on Fig. [U The 
integral /* (see caption for a and b) was calculated by using Gauss-Legendre quadrature with 10, 20, 
50, 100, 200, 500, 1 000, 2 000, and 5 000 points. Starting from the second step {k - 2, i.e. 20 points), 
the relative discrepancy was calculated, 

|Re(4%-Re(4Vi| ^ 

|Re(4\l ^ |Im(4'),| 

where (7*)^ is the value of /*, obtained on the ^-th step. For the sake of illustration, here we were 
using simple Yamaguchi-style real form factor with constant 2.7 fm '. It should be noted, that for 
this large value of 77 = 5.3 and small value of o' = 0.5 fm“* (the corresponding system is p + 
at Ec.m. - 6 MeV) the final result is being calculated by subtraction of two large numbers, and the 
integral 7* must be computed with crt ^ 10 to get even the first digit of the final result. This new 
regularization scheme will be essential for successful progress in the numerical implementation of the 
Faddeev-AGS equations in Coulomb basis including spin degrees of freedom. 
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